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Abstract 

In this article, we calculate the contribution of the higher-twist Feynman diagrams to the large- 
Pt inclusive single pion production cross section in photon-photon collisions in case of the running 



coupling and frozen coupling approaches within holographic QCD. We compare the resummed 



higher-twist cross sections with the ones obtained in the framework of the frozen coupling approach 
and leading-twist cross section. Also, we show that in the context of frozen coupling approach a 
higher-twist contribution to the photon-photon collisions cross section is normalized in terms of 
the pion electromagnetic form factor. 
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I. INTRODUCTION 



One of the most significant theoretical advances in recent years has been the apphcation of 
the AdS / CFT correspondence jll between string theories defined in five- dimensional anti- 
de Sitter (AdS) space-time and conformal field theories in physical space-time. Quantum 
chromodynamics (QCD) is not itself a conformal theory; however there are indications, both 
from theory 2, 3] and phenomenology 0,0] that the QCD coupling is slowly varying at small 
momentum transfer. In addition, one can argue that if the gluon has a maximum wavelength 
or an effective mass due to confinement, that gluonic vacuum polarization corrections 
and the /3-function must vanish in the infrared. If there is a conformal window where the 
QCD coupling is large and approximately constant and quark masses can be neglected, 
then QCD resembles a conformal theory, thus motivating the application of AdS/QCD to 
QCD. So, even though QCD is not conformally invariant, one can use the mathematically 
representation of the conformal group in five- dimensional anti-de Sitter space to construct 
an analytic first approximation to the theory. The resulting AdS/QCD model gives accurate 
predictions for hadron spectroscopy and a description of the quark structure of mesons and 
baryons which has scale invariance and dimensional counting at short distances, together 
with color confinement at large distances. 

The hadronic wave function in terms of quark and gluon degrees of freedoms plays an 
important role in QCD process predictions. For example, knowledge of the wave function 
allows to calculate distribution amplitudes and structure functions or conversely these pro- 
cesses can give phenomenological restrictions on the wave functions. 

In References jzHl the higher- twist effects was calculated within the frozen coupling 



constant approach. In Ref. [12!], it was noted that in perturbative QCD (pQCD) calculations, 
the argument of the running coupling constant in both renormalization and factorization 
scale should be taken as equal to the square of the momentum transfer of a hard gluon 
in a corresponding Feynman diagram. But defining in this way, as(Q^) suffers from infrared 
singularities. 

The contribution of large orders of perturbation theory related to the socalled renor- 



malons has 
summation 



jeen investigated by several authors, using in particular the method of Borel 
13l-ll7l|. In the case of (QCD), the coefficients of perturbative expansions in 



the QCD coupling can increase dramatically even at low orders. This fact together with 



2 



the apparent freedom in the choice of renormahzation scheme and renormahzation scales, 
hmits the predictive power of perturbative calculations, even in applications involving large 
momentum transfer, where is effectively small. 

Investigation of the infrared renormalon effects in various inclusive and exclusive pro- 
cesses is one of the most important and interesting problems in the perturbative QCD. 
It is known that infrared renormalons are responsible for factorial growth of coefficients 
in perturbative series for the physical quantities. But, these divergent series can be re- 
summed by means of the Borel transformation 13] and the principal value prescription [l^ 
and effects of infrared renormalons can be taken into account by a scale-setting procedure 
OisiQ"^) '^s(exp(/((5^))<5^) at the one-loop order results. Technically, all-order resum- 
mation of infrared renormalons corresponds to the calculation of the one-loop Feynman 
diagrams with the running coupling constant A;^) at the vertices identically equivalent 
to calculation of the same diagrams with nonzero gluon mass. 



In this work, we apply the running coupling approach 19| in order to compute the 
effects of the infrared renormalons on the pion production in photon-photon collisions within 



holographic QCD. This approach was also employed previously |20l424i| to calculate the 
inclusive meson production in photon-proton, proton-proton and photon-photon collisions. 

Since experiments examining high-p^ particle production in two-photon collisions have 
been improved, it becomes important to reassess the various contributions which arise in 
quantum chromodynamics. Also the experimental measurement of the inclusive charged 
pion production cross section is important for the photon-photon collisions program at the 
future International Linear Collider (ILC). 

Therefore, the calculation and analysis of the higher-twist effects on the dependence of 
the pion wave function in single pion production at photon-photon collisions by the running 
coupling approach within holographic QCD are very interesting search points. 

In this respect, the contribution of the higher-twist Feynman diagrams to a single meson 
production cross section in photon-photon collisions is computed by using various pion wave 
functions from holographic QCD. Also, the leading and resummed higher-twist contributions 
are estimated and compared to each other. 

We organize the paper as the follows; In Section [Tll we provide some formulae for the 
calculation of the contributions of the higher twist and leading twist diagrams. In Section 
\\\\\ we present some formulae and analysis of the higher twist effects on the dependence of 
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the pion wave function by the running couphng constant approach, and in Section llVt the 
numerical results for the cross section and discussion for the dependence of the cross section 
on the pion wave functions are presented. Finally, some concluding remark are stated in 
Section |Vl 



II. HIGHER TWIST AND LEADING TWIST CONTRIBUTIONS TO INCLU- 
SIVE REACTIONS 

The higher-twist Feynman diagrams for the pion production in the photon-photon colli- 



sion 77 — )■ MX are shown in Fig. 1(a). T 



means of the Brodsky-Lepage formula 



le amplitude for this subprocess can be found by 



25| 



M(s,t)= / dXi / dx25{l - Xi - X2)^m{,Xi,X2jQ )TH{s,t]Xi,X2). (2.1) 

Jo Jo 

In Eq.(2.1), Th is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. The hard-scattering amplitude Th-(s, t; xi, X2) depends on a process and can 
be obtained in the framework of pQCD and represented as a series in the QCD running 
coupling constant ^^(Q^). The light-cone momentum fractions = 1 — X specify 

the fractional momenta carried by quark and antiquark in the Fock state. As higher-twist 
subprocess which contribute to 77 — )■ vrX, we take ■yq — )■ Mq. 

The Mandelstam invariant variables for subprocesses 'jq — )■ Mq are defined as 

s = {pi+ P-if, i = (p^ - pmY, u = {pi- pmY- (2.2) 

We have aimed to calculate the pion production cross section and to fix the differences due 



to the use of various pion model wave functions. The asymptotic pion wave functions |26l.l27| 
and the Vega-Schmidt-Branz-Gutsche-Lyubovitskij (VSBGL) [28]wave function predicted by 
AdS/QCD, and also the pQCD evolution has the form: 



^SBcd^) = ^vW^)exp I ) , <l>Lyi^) = V3Ux{^-x) (2.3) 

where = 92AMeV is the pion decay constant. 

We now incorporate the higher-twist subprocess 7g — )■ Mq into the full inclusive cross 
section. In this subprocess photon and the pion may be viewed as an effective current 



striking the incoming quark line. Therefore, the complete cross section in formal analogy 
with deep-inelastic scattering, is written as 



diCT 3 d(j 

E-^{-f-f ^ MX) ^ dx5{s + i+u)sGg/^{x,-i)^{'yq^ Mq) + (t^u). (2.4) 

m 

Here Cqj^ is the per color distribution function for a quark in a photon. The subprocess 
cross section for tt, and pr production is 

da f ^^^^[L>(s,^i)]2 1 [4 + 4] , M = 7r,pi, 



where 



D{s,u) = eiil^{Ql)as{Ql) + e2uh{Ql)as{Ql), 



(2.6) 



dx 



Q^) 

x{l — x) 



(2.7) 



and Q\ — s/2, Q\ — —u/2, represents the momentum squared carried by the hard gluon 
in Fig. 1(a), 61(62) is the charge of ^1(^2) ^-nd Cp — 

The factors reflect the exclusive form factor of the pion as is the motivation the argu- 
ments of as and /jr. Note that the relation between and the pion form factor completely 
fixes the normalization of the higher-twist subprocess. The full cross section for tt and 
production is given by 



da , , s 

d'^p s -\-u 



^Gg/^{x,-i) 



3 



99 



P ~'~ u'^ 



+ 



SnaECp [D{s,i)f 
S + t^ .//V ' /' § s^(-u) 

qq ^ ' 



1 1 



(2.8) 



As seen from Eq.(2.8), the subprocess cross section for longitudinal pl production is very 
similar to that for tt production. We have extracted the following higher-twist subprocesses 
contributing to the two covariant cross sections in Eq.(2.5) as 



791 (91^2)52 , 752 (9192)52- 



(2.9) 



Also from Eq.(2.8), at fixed pr, the cross section falls very slowly with s. Additionally, 
at fixed s, the cross section decreases as multiplied by a slowly varying logarithmic 



function which vanishes at the phase-spase boundary. Thus, the pr spectrum is fairly inde- 
pendent of s except near the kinematic hmit. 

One of the important problem in the single inclusive pion production in photon-photon 
collision is the possibility of normalization of the higher-twist subprocess cross section in 
terms of the electromagnetic form factor of the pion. The electromagnetic form 

factor Ft^{Q'^) of pion is given by the expression 

FAQ^)= [ [ dxidx2K{xi,X2,Q^)TH{xi,X2,as{XQ^),Q^)<^>^{x,,X2,Q^). (2.10) 
Jo Jo 

This allows us to completely determine the ^yq — > Mq cross section in terms of the pion form 
factor, through the relation 

iim (2.11) 



47rC,.a,(g2) 

It should be noted from Eq.(2.11) that the form factor contains the square of /^(Q^)- In 
principle, experimental measurement of F^((5^) determines and hence the 7g — )• Mq. So 
one can determine the cross section of 77 — > MX explicitly. 

Now wc can conclude that in the frozen coupling constant approach vr meson production 
higher-twist cross section of 77 MX is normalized in terms of the pion electromagnetic 
form factor. 

Extracting the higher-twist corrections to the pion production cross section and a com- 
parison of higher-twist corrections with leading-twist contributions are essential problems. 
The contribution from the leading- twist subprocess 77 — > qq is shown in Fig. 1(b). The 
corresponding inclusive cross section for production of a meson M is given by 
da 



d^p 



- E / ^^(^ + i + ^ (2-12) 

- Jo z at 



where 

t ^ u 
s — s, t — - u — —. 



z z 



Here s, t, and u refer to the overall 77 MX reaction. Dg {z, —i) represents the quark 
fragmentation function into a pion containing a quark of the same flavor. For vr^ produc- 
tion we assume that Dt,+ = In the leading-twist subprocess, pion is indirectly 
emitted from the quark with fractional momentum z. The final form for the leading- twist 
contribution to the large-pT pion production cross section in the process 77 — > MX is 



u i 



(2.13) 



The contributions from these leading-twist subprocesses strongly depend on some phe- 
nomenological factors, such as, quark and gluon distribution functions in meson and frag- 
mentation functions of various constituents etc. We should note that D{z, —i)/ z behaves 
as l/z"^ when ^ ^ 0. For the kinematic range considered in our numerical calculations, 
D{z, —i)/z increases even more rapidly. We obtain of the final cross section, Eq.(2.13), as 
follows: At fixed p^, the cross section decreases with s asymptotically as 1/s. At fixed s, the 
D{z, —t) function causes the cross section to decrease rapidly when px increases towards the 
phase-spase boundary {z ^ 1). As s increases, the phase-spase boundary moves to higher 
Pt, and the pr distribution broadens. 



III. HIGHER TWIST EFFECTS WITHIN HOLOGRAPHIC QCD AND THE 
ROLE INFRARED RENORMALONS 

The main problem in our investigation is the calculation of the integral in Eq.(2.6) by the 
running coupling constant approach within holographic QCD. It should be noted that, in the 
exclusive processes, the coupling constant runs not only due to loop integration, but also 
because of the integration in the process amplitude over the light-cone momentum fraction of 
hadron constituents. In this respect, the exclusive processes have two independent sources 
of power corrections to their characteristics. One of the loop integration and the second 
integration over the light-cone momentum fraction of hadron constituents. Therefore, it is 
worth noting that, the renormalization scale (argument of ag) according to Fig. 1(a) should 
be chosen to be equal to = {1 — x)s, and Q2 = —xu. The integral in Eq.(2.6) takes the 
form in the framework of the running coupling approach 

The as(A/i|jj,) has the infrared singularity at x — >^ 1, for A 1 — x or x — > 0, for A = x and 
so the integral (3.1) diverges. For the rcgularization of the integral, we express the running 
coupling at scaling variable q;s(A;U^^^) with the aid of the renormalization group equation 
in terms of the fixed one as{Q'^)- The solution of renormalization group equation for the 
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running coupling a = cts/vr has the form 

a(A) 



a 



1 + In A 
4 



Then, for as{\Q^), we get 



a. 



l + \n\/t 
where t = 47r/a^((5^)/3o = 4/a/3o. 

Having inserted Eq.(3.3) into Eq.(2.6) we obtain 

\ a,{Xfil)<i>M{x,Ql] 



(3.2) 



(3.3) 



D(s,u) = eiu I dx^^ — ^ — \ ' + 628 I dx 



eiuas{s)ti / dx 



x{l-x) 
<l>Mix,Ql) 



\ as(A/i|J$M(a;,Qi) 



+ e2sas{—u)t2 I dx 



x{l — x) 



x(l - x)(t2 + In A) 



(3.4) 



/g x(l — x){ti + In A) 
where ti = 4:n/as{s)f3o and ^2 = 47r/a;s(— 'u)/3o- 

Although the integral (3.4) is still divergent, it is recast into a suitable form for calculation. 
Making the change of variable as z = In A, we obtain 



.X . .... <i>M{x,Q^)dx , „ /"i $M(x,g2)c?a; 

19(5, m) = eiMtts^sjti / — — ; — - + e2sas[—ujt2 



Q x{l-x){ti + z) 



X{l-X){t2 + Z) 



(3.5) 



In order to calculate (3.5), we will apply the integral representation of l/(t + z) [30|, |3l| as 



t + z 



(3.6) 



which gives 



iJ(s, -u) = eitos(s)ti / / ^ \-e2uas{—u)t2 



^0 



x(l — x) 



~'0 



x(l — x) 
(3.7) 



In the case ^afyi^) for the D{s,u), Eq.(3.7) is written as 



D(s,u) 



"'0 



2 2 







due-'--B(^-^--v^ (3.8) 



and for ^^^^(a;) wave function 



4y37r/, 



eiu 



f3o 



due 



-tiu 



1-U 



+ 



628 



/3o 



due 



-t2U 



l-u 



(3.9) 



where B{a,(3) is Beta function. The structure of the infrared renormalon poles in Eq.(3.8) 
and Eq.(3.9) strongly depend on the wave functions of the pion. To remove them from 
Eq.(3.8) and Eq.(3.9) we adopt the principal value prescription. We denote the higher-twist 
cross section obtained using the running coupling constant approach by (E^^)'^'^*. 
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IV. NUMERICAL RESULTS AND DISCUSSION 



In this section, we discuss the numerical results for higher-twist and renormalon effects 
with higher-twist contributions calculated in the context of the running coupling constant 
and frozen coupling approaches on the dependence of the chosen pion wave functions in 
the process 77 — MX within holographic QCD. For the higher-twist subprocess, we take 
7*?! ~^ (91^2 )Q'2) and 7^2 — ^ (q'i^2)Q'2 contributing to 77 — )■ MX cross sections. Inclusive pion 
photoproduction represents a significant test case in which higher-twist terms dominate 
those of leading-twist terms in certain kinematic domains. For the dominant leading-twist 
subprocess for the meson production, we take the photon-photon annihilation 77 — )■ qq in 
which the vr pion is indirectl y e mitted from the quark. The quark distribution function 
inside the photon was used [32| and the gluon and quark fragmentation functions into a 



pion was used [33 1. 

The results of our numerical calculations are plotted in Figs. 2-15. Firstly, it is very 
interesting comparing the higher twist cross sections obtained within holographic QCD with 
the ones obtained within perturbative QCD. In Fig. 2 and Fig.3 we show the dependence of 
higher-twist cross sections {T^^'^Y , and {Ti^+Y^'^ calculated in the context of the frozen and 
running coupling constant approaches as a function of the pion transverse momentum px for 
<l>^°'(x), and ^v'sbgl{^) pion wave functions at ?/ = 0. It is seen from the figures that 

the higher-twist cross section is monotonically decreasing with an increase in the transverse 
momentum of the pion. In Fig.4-Fig.7, we show the dependence of the ratios {11^+ Y, 

{T.^^Y^'/iT.^^f, (Ef7)7(S^^) and (EfJ)^'=V(E^J) as a function of the pion transverse 
momentum px for ^^{x) and ^^sBGhi^) pion wave functions. Here is the 

leading-twist cross section. As shown in Fig. 4, in the region 20 GeV/c < < 80 GeV/c 
higher- twist cross section for is suppress by about half orders of magnitude relative to 

the higher-twist cross section for ^v'sbgl^^)^ regions 10 GeV/c < pt < 20 GeV/c 

and 80 GeV/c < px < 90 GeV/c, higher-twist cross section is suppress by about two 

orders of magnitude relative to the higher-twist cross section for {^v'sbgl)^- Also higher- 
twist cross section for is suppress by about half orders of magnitude relative to the 
higher-twist cross section for ^v°SBGLi^) ■ Fig-5, Fig.6 and Fig. 7, the dependence of the 
ratios (Sf^)^^^V(SfJ)0, (SfJ)7(S^J) and (Sf7)''^7(S^^) are displyed as a function of the 
pion transverse momentum p^ for the and ^v'sbglI^) pion wave functions. It 
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is see that, resummed higher-twist cross section for is suppress by about two orders 

of magnitude relative to the higher-twist cross section for {^iFsbglT- Noticed that one- 
half order is suppress for (E^°^)'' and one order is suppress for (E^^j^)°. It is observed from 
Fig.6 and Fig.7 that, the ratios (SfJ)7(S^J) and (SfJ)''"7(S^;f^) for all wave functions 
decrease with an increase in the transverse momentum of pion. In Fig.8 - Fig. 10, we have 
depicted higher-twist cross sections {^'i^^f, and ratio (E^";^)7(Ef7)°, (E^°|)''^7(EfJ)°, as 
a function of the rapidity y of the pion at 7^ — 183 GeV and px — 14.6 GeV/c. The figures 
show that the higher-twist cross section and ratios have a different distinctive behavior. 
As is seen in Fig.9 ratio i^afy)^ / (^v'sbglT has a maximum approximately at the point 
y = —1.92. However in this point ratio {'^v'sbgl)^ /i^asy)^ ^ minimum. As is seen 
from Fig. 10 resummed higher-twist cross section for is suppress by about one order 

of magnitude relative to the resummed higher- twist cross section for ^asyi^) with an 
increasing rapidity of pion ratio is kepped approximately constant. But resummed higher- 
twist cross section for ^^°y(x) is suppress by about one half- two orders of magnitude relative 
to the frozen higher-twist cross section for ^v'sBGii^) ^ maximum approximately at 

the point y = —1.92. However resummed higher- twist cross section for is suppress 

by about one order of magnitude relative for resummed higher-twist cross section for 
and to stay is constant with an increasing rapidity of pion. Figures also show that, the ratio 
depends on the choice of the pion wave function. Analysis of our calculations concludes that 
(E^_^)°, and {T,^+Y^^ higher-twist cross sections and ratio sensitive to pion wave functions 
predicted by holographic and perturbative QCD. 

We have also carried out comparative calculations in the center-of-mass energy 7^ = 
209 GeV and obtained results are displayed in Fig.ll-Fig.15. Analysis of our calculations at 
the center-of-mass energies ^/s = 183 GeV and ^/s — 209 GeV, show that with increasing in 
the beam energy contributions of higher twist effects to the cross section decrease by about 
1-2 orders of magnitude. As is seen from Fig.4, Fig. 5, Fig.9, Fig. 10, Fig.l2, Fig. 14 and Fig. 15 
that infrared renormalon effects enhance the perturbative predictions for the pion production 
cross section in the photon-photon collisions about 1-2 orders of magnitude. Our opinion 
is this feature of infrared renormalons may help the explain theoretical interpretations with 
future experimental data for the pion production cross section in the photon-photon colli- 
sions. In our calculations of the higher-twist cross section of the process the dependence of 
the transverse momentum of pion appears in the range of {10~^-^10~'^'^)mb/GeV'^. Therefore, 
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higher-twist cross section obtained in our work should be observable at ILC. 

V. CONCLUSIONS 

In this work, the single meson inclusive production via higher twist mechanism within 
holographic QCD are calculated . For calculation of the cross section the running coupling 
constant approach is applied and infrared renormalon poles in the cross section expres- 
sion are revealed. Infrared renormalon induced divergences is regularized by means of the 
principal value prescripton and the Borel sum for the higher twist cross section is find. It 
is observed that, the resummed higher- twist cross section differs from that found using the 
frozen coupling approximation in some region considerably. We proved that the higher-twist 
cross section for n pion production in the photon-photon collisions may be normalized in 
terms of the pion form factor. The following results can be concluded from the experiments; 
the higher-twist contributions to single meson production cross section in the photon-photon 
collisions have important phenomenological consequences, the higher-twist pion production 
cross section in the photon-photon collisions depends on the form of the pion model wave 
functions and may be used for their study. Also that the contributions of renormalons effects 
within holograpich QCD in this process is essential and may help to analyse experimental 
results. Further investigations are needed in order to clarify the role of higher-twist effects 
in QCD. Especially, the future ILC measurements will provide further tests of the dynamics 
of large-pT hadron production beyond the leading twist. 
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FIG. 1: (a): The higher-twist contribution to 77 — )• MX] (b): The leading-twist contribution to 

77 MX 
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FIG. 2: Higher-twist vr"'" production cross section (Ti^^)^ as a function of the pT transverse mo- 
mentum of the pion at the cm. energy ^/s = 183 GeV. 
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FIG. 3: Higher-twist vr^ production cross section (S^^^)'"'^^ as a function of the px transverse 
momentum of the pion at the cm. energy ^/s = 183 GeV. 
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FIG. 4: Ratio (S^|,)^/(S^J')'', where higher-twist contribution are calculated for the pion rapidity 
y = at the cm. energy ^/s = 183 GeV as a function of the pion transverse momentum, pj-- 




FIG. 5: Ratio {H'^^^Y'^''^ / {^^^)^ , as a function of the pt transverse momentum of the pion at the 
cm. energy ^/s = 183 GeV. 
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FIG. 6: Ratio (S^)^/(E^+), where higher-twist contribution are calculated for the pion rapidity 
y = at the cm. energy ^/s = 183 GeV as a function of the pion transverse momentum, pj-- 
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FIG. 7: Ratio {T,^+y'^''^/{^^+), as a function of the pt transverse momentum of the pion at the 
cm. energy ^/s = 183 GeV. 
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FIG. 8: Higher-twist vr"'" production cross section (S^j,)^ , as a function of the y rapidity of the 
pion at the transverse momentum of the pionpr = 14.6 GeV/c, at the cm. energy ^/s = 183 GeV. 
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FIG. 9: Ratio (S^j.)'^/(S^4^)'', as a function of the y rapidity of the pion at the transverse 
momentum of the pion px = 14.6 GeV/c, at the cm. energy y/s = 183 GeV. 
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FIG. 10: Ratio as a function of the y rapidity of the pion at the transverse 

momentum of the pion px = 14.6 GeV/c, at the cm. energy y/s = 183 GeV. 
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FIG. 11: Higher-twist vr"'" production cross section as a function of the pT transverse 

momentum of the pion at the cm. energy y/s = 209 GeV. 
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FIG. 12: Ratio {'E'^Jj^Y'^^ / {T,^+)^ , as a function of the px transverse momentum of the pion at the 
cm. energy ^/s = 209 GeV. 
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FIG. 13: Higher-twist vr"*" production cross section (S^j,)'^ , as a function of the y rapidity of the 
pion at the transverse momentum of the pion px = 16.7 GeV/c, at the cm. energy t/s = 209 GeV. 
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FIG. 14: Ratio as a function of the y rapidity of the pion at the transverse 

momentum of the pion px = 16.7 GeV/c, at the cm. energy y/s = 209 GeV . 
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FIG. 15: Ratio {Ti^+Y^'^/{Ti^^)^, as a function of the y rapidity of the pion at the transverse 
momentum of the pion px = 16.7 GeV/c, at the cm. energy y/s = 209 GeV . 
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